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Abstract. Three types of integral representations for the cumulative distribution 
functions of convolutions of T p (oik, Efc, A^)-distributions with non-centrality matri- 
ces Afe are given by integration of products of simple complex functions over the 
p-cube (— 7r,7r] p . In particular, the joint distribution of the diagonal elements of a 
generalized quadratic form XAX' with n independent A/j,(/ifc, E)-distributed columns 
in X p xn and a fixed A > is obtained. For a single F p (a, £, A)-cdf (p — l)-variatc 
integrals over (— 7r,7r] p_1 are derived. The integrals are numerically more favourable 
than integrals obtained from the Fourier- or Laplace inversion formula. 
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1. Introduction 

The following notations are used: E(n) stands for ^ ni+ + n —„ with 
rii, . . . , rip G No and ^ without any indices means X)^LoS(n)- The notation D > 
is also used for non-symmetrical matrices D pxp with only non-negative eigenvalues. 
The spectral norm of a p x p-matrix B is denoted by ||S|| , / or I p is always an identity 
matrix and C p is the p-cube (— ir, tt] p . 

The Laplace transform (L.t.) of a p-variate non-central T p (a, S, A)-density with 
a > 0, £ > and a non-centrality matrix A > was originally obtained from the L.t. 
of a non-central W p (2a, E, A)-Wishart distribution (with an additional scale factor 2) 
and is given by 

/(ir,..., t p ; a, £, A) = \I p + ET|- Q etr(-ET(7 + ST)" 1 A), (1) 
T = diag^i,...,^),*!,...,^ > 0. 

This function / is generally the L.t. of the density of a real measure on (0,oo) p which 
is not always a probability measure. The term "T p (a, E, A)-distributkm" is used here 
in this general sense. The exact set of values a, leading to a probability density (pdf) 
f(x±, . . . , x p ; a, E, A), depends on E and presumedly on A. To obtain a pdf, all positive 
integers 2a (degrees of freedom) are admissible and all 2a > p — 1. Moreover, in the 
central case all non-integer values 2a > p — 2 > are allowed. For p — 2<2a<p— 1 
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see Royen (1997). Furthermore all a > are admissible if \I + ST| _1 is infinitely 
divisible. Two characterizations of infinite divisibility of a T p (a, E)-distribution are 
found in Griffiths (1984) and Bapat (1989). Further conditions for admissible non- 
integer 2a < p — 2 are given in Royen (1997), (2006). 

Three integral representations by integration over C p are provided by theorem 2 in 
section 4 for the functions 

F(x!, . . . ,ar p ; ai, . . . , a n , Si, . . . , £ n , Ai, . . . , A„) (2) 

= I •••/ • ■ • ,£ P ;ai) ■ • ■ ,a„,Si, . . . ,E„, Ai, . . . , A„)d£i . . .d£ p , 

Jo Jo 

where / has the L.t. 

n 

Y[ \I P + E fc Tr Qfc ctr(-£ fc T(/ + SfcTJ-^fc), (3) 
fe=i 

ai, . . . , a„ > 0, Ei, . . . , E„ > 0, Ai, . . . , A„ > 0. 

Thus, F is not always the cumulative distribution function (cdf ) of a probability measure. 

In particular let X pxn be &Af pxn (M pxn , Ep Xp (g)/ n ) -random matrix and A nxn > of 
rank q with T' AT = A = diag(Ai, . . . , A„), Ai > . . . > A„. Then the joint distribution of 
the diagonal elements of the generalized quadratic form \XAX' equals the distribution 
of the diagonal of \YKY' with &Af pxn (MT, E pxp ®/„)-distributed Y = XT. This is the 
distribution of a sum of q independent r p (|,AfeE,Afc = \^%^* k E _1 )-random vectors, 
where \i\ is the fc-th column of M* = MT. This joint distribution of p quadratic 
forms of normal random vectors is comprised within theorem 2 as a special case with 
ct/c = i, E fe = A fe E, k = 1, . . . , q. For methods under more general assumptions see also 
Blacher (2003). For a survey of univariate quadratic forms of normal random variables 
see chapter 4 in Mathai and Provost (1992). For several quadratic forms of skew elliptical 
distributions see B.Q. Fang (2005). 

In Royen (1991), (1992) three different types of series expansions for the Xp(2a, E)- 
cdf were derived from three different representations of the X P (2a, E)-L.t. which are 
extended to the general T p (a, E, A)-L.t. in section 3 in a similar way as in Royen 
(1995). 

Some series expansions, closely related to the first two types, are already found in 
Khatri, Krishnaiah and Sen (1977). The third type was introduced because of its superior 
convergence properties. The simple method to transform many series expansions into 
integrals over C p is explained in more detail in section 2 and summarized in theorem 1. 
The idea is as follows: 

If A(z\, . . . , z p ) and B{z\ 1 . . . , z p ) are analytical functions whose power series have 
the coefficients a (mi, ... , m p ) and b(n\, . . . , n p ) and which are absolutely convergent for 
max|2j| < ta and maxlz^l < tb respectively, where r^ 1 < ta, then 

(2tt)-p / A(y 1 ,...,y p )B(y^ 1 ,...,y- 1 )dip 1 ...di Pp (4) 

= ^2 a ( ni ' ■ • ■ ' n p) & ("i> ■ • ■ ' n p) 

holds with yj — re lip i , — 7r < (fj < ir, j = 1, . . . ,p and r^ 1 < r < Ta- 
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The integrals in J?]) might be more economical than the series if the generating 
functions A and B are simple available functions and if the series are slowly convergent 
with very intricate coefficients. For non-central multivariate gamma distributions series 
expansions are practically not feasible. 

The integral representations in theorem 2 of section 4 are of the type in (|4|). As 
long as no elementary density formulas are availale it should be a reasonable way to 
obtain the joint cdf by integration of elementary terms only over C p and not over W as 
by the Fourier or Laplace inversion formula. A single T p (a, E, A)-cdf is represented by 
a (p — l)-variate integral over C p _i in section 5. 

A totally different ( m ^~ )— variate integral representation of the T p (a, E)-cdf has 
been given recently by Royen (2006), which is based on m-factorial decompositions 
^2 p x p = D — BB' , where D is a real or complex diagonal matrix minimizing the rank 
m of E -1 — D. Approximations to a T p (a, E)-cdf are obtained by m-factorial approx- 
imations to E with a low value of m. These approximations are improved further by 
successive correction terms. 



2. The method 



Theorem 1 in this section can be generalized in many ways, e.g. for Fourier trans- 
forms, but the version below is sufficient for the purpose of the underlying paper. 

Let f(t\, . . . , t p ), t±, . . . , t p > 0, be a given L.t. of an unknown function f(x\, . . . , x p ) 
with / = for miniEj < 0. It is assumed that there are univariate L.t. 9j {t) of 
some probability densities gj (x) on (0, oo) and further functions hj(t) with < 1, 

uniformly for t > 0, which enable a representation 

/(*!,..., t p ) = {^^ j0 {t 3 ^B{h x {h),...,h p {t P V (5) 

with an analytical function B{z\, . . . , z p ) whose power series expansion 

p 

^ & ( n i>-- Y[ z T ( 6 ) 

3=1 

is absolutely convergent for \z±\, . . . , \z p \ < tb with a certain value rs > 1. 

Furthermore, the products gj (t)(hj(t)) n are supposed to be the L.t. of continuous 
functions gj n (x), x > 0, which satisfy the conditions 

< n c k(x) with a constant c and 

/•oo {') 

/ k(x)e~ tx dx < oo for all £ > . 
Jo 

Hence, the generating functions (generators) 

oo 

9j(x,y) = J £2g j n(x)y n , j = l,...,p, (8) 
are defined for all x > and \y\ < 1, and they have the L.t. 
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Theorem 1. Under the assumptions from and / in |5|) is the L.t. of 

r p 

f(x lt ...,x p ) = (2n)-* Bfc\...,y; 1 )]lg J (x j ,y j )d<p j (10) 

J C P j=l 

with yj = re lVi , —it < <p <tt, r^ 1 < r < 1, gj from f5|). 
Proof. The integral in (fTU|) is evaluated by 

( 27r r p J c (^2 K m U ■ ■ ■ , Tip) JJ y7™ J j J| ffjn, Oj)y" 3 ' j d^i . . . 

= J2 b{nx, ■ ■ ■ ,n p ) n^=i ftnj (acj) 
and this series has the L.t. from ([5]). 
Some further remarks: With 

G j( x j>yj)= 9j(€,yj)d£ (ii) 

Jo 

instead of the gj in dTDJ) , the corresponding representation arises for 

F( Xl ,...,x p )= H ... H /& (12) 

Jo Jo 

If the series in ([5]) are absolutely convergent for all y £ C then additionally 

Jim hj(t) = (13) 



t — >oo 



is supposed to hold. Then the rhs of §§§ is the L.t. of gj(x,y) for any fixed y and all 
sufficiently large t. 

In some cases the functions gj and their L.t. cjj are known from univariate marginal 
distributions apart from some scale factors. If the functions Uj — hj(t) are explicitly 
invertible then 

B( Ul ,...,u p ) = — - (14) 

can sometimes be found easily from the given /. 

3. Three representations for the T p (a, E, A)-Laplace transform 
and the related generators 

With any v > we define 

Zj = (1 + v~ tj)~ , <j > 0, ttj = 1 — Zj = v~ 1 tjZj, uij = zj — Uj, 
Z = diag(zi, . . . , z p ), U = diag(ui, . . . ,u p ), f2 = diag(wi, . . . , u p ). 



(15) 
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The scale factor v is introduced to obtain ||_B|| < 1 for the matrices B defined in (f5U|) 
below and to effect the convergence of some series expansions. For a more general scaling 
see remarks following theorem 2 in section 4. 
From the relations 

v~ l T = UZ-\ l p = Z + U, Q = Z-U, (16) 

it follows for the matrices / + ST in the L.t. fTJ): 

I + Y<T = I + VY.UZ- 1 = {Z + vY,U)Z- 1 (17) 

and 



I+(vH- I)U, (18a) 
Z + vEU = { vE(I + (ir 1 ^ 1 - I)Z), (186) 
i(7 + uS)(J+ (2(1 + (18c) 



and therefore 



ST|- Q = c a \Z\ a \I + BY\- a (19) 



with 



Y = U, B = vT,-I, c = 1, (20a) 

Y = Z, B = (vY)- 1 -I, c=\I + B\, (206) 

Y = fl, B = 2(1 + VT,)- 1 -I, c=\I + B\. (20c) 

It should be noticed that ||B|| < 1 in (j2"0b ) for every v > and E > 0. 

Now, using (fT6|) , by a straightforward calculation the L.t. in can be represented 

by 

f(t 1 ,...,t p ;a,T,,A) = 



( \Z\ a \I + BU\- a etr(-(I + B)U(I + BU^A), (21a) 

I \I + B\ a eti(-A)\Z\ a \I + BZ\- a etr(Z(I + BZ)- 1 ^ + B)A), (21b) 

\I + B\ a etr(-±A(I - B)) (21c) 
■\Z\ a \I + SO|- a etr(|fi(7 + Bn)- x (I + B)A(I - B)), 

with the corresponding matrices B from (|20|) and Z, U, fl from ([15]). 
For the former series expansions the following relations were used: 
Laplace transform f(t): f( x )'- F{ x ) — Jo f(0^ : 



vg£> n (vx) C( ra ) ( 2 2a) 

z a+n vg a+n (vx) G a+n (vx) (226) 

z a uj n vh a , n (vx) H a>n (vx) (22c) 



where z = (1 + y-H)' 1 , g a +n{x) = e' x x a ' 1+n /T(a + n), 

L^\x)g a {x) 



(n) / \ d n (a-l + n 



with the generalized Laguerre polynomials L„ 1 and 



= (-1)" 



a — 1 + n 
n 



L^ 1 \2x)g a {x). 



The last identity is verified by L.t. 

The following bounds are derived from (22.14.13) in Abramowitz and Stegun (1965): 



(™) / \ 

9a+n\ x ) 



< 



e x ^ 2 g a (x), a > 1 

2na~ 1 e x / 2 g a (x), < a < 1 



(23) 



1 ' Wl " | 2nx a - 1 /T{a + l), 0<a<lj' y ' 

matching with the conditions in ([7]). 

The following generators (generating functions) with the r(a + n)-cdf G a + n (x) are 
required for the formulas in theorem 2: 

[ EZo G { :Ux)y- = j^G a (x, , \y\ < 1, (25a) 
F a (x,y) = I Y.n=vG a +n{x)y n = G a (x,y), y e C, (256) 

[ EZo H a , n {x)y n = j^G a (x, ^) , \y\ < 1 (25c) 

The identities (a) and (c) are verified by the L.t. of f a (x,y) = -j^F a {x,y). A short 
calculation shows 

' (G a (x) - y^e^y- 1 ^ G a (xy)) , y ? 1, a>0 

= (Ga-i(^) - y 1 -^^- 1 ) 1 G^!^)) , a>l, G := 1 (26) 

^ xg a (x) + (1 +x - a) (^(a;), y = l 

and 

9 fgaW + y 1 -^- 1 )* G„(ay), «>0 1 

^( a; ^) = fe G «^y) = ( yl - Qe (,- 1 ), GQ _ l(xy); a >! )■ 

The functions F a (x, y) are especially simple for a e N since G Q (z) = 1 — S?=o zJ /j'; 
a S N. 

Besides, 

G k+1/2 (z) = erfiz 1 / 2 ) - e~'J2 * G N °' 

The following simple lemma is used for the proof of theorem 2. 
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Lemma 1. If B is a symmetrical p x p-matrix with \\B\\ < 1 and 
Y = diag(?/i, . . . , y p ) then the power series expansion 

p 

\I + BY\- a =J2b(ni,...,n p )H 



3 = 1 



is absolutely convergent for max|?/, | < re = \\B\ 



This follows from Y2( n ) IK 71 !) ■ ■ ■ i n p)\ — 0(-d n ) with any d > \\B\\, which has been 
already shown in (2.1.16) . . . (2.1.18) in Royen (1991) (with the notation — C instead of 
B). 

4. The integral representations 

In theorem 2 below the functions F{x\, . . . , x p ; a\, . . . , a n , Si, ... , £ n , Ai, . . . , A„) 
from @ are represented by three different integrals over C p = (— tt,tt] p . Together with 
the generators F a from (|25[) . a = Y], —1 a k , the following matrices are used with a scale 
factor v to enforce \\B k \\ < 1: 

B k = vX k - I, D k = A k (I + B k ), F a from (27a) 

B k = {v^ k )- 1 -I, D k = (I + B k )A kl F a from (2Sb), (276) 

B fc = 2(/ + «S fc )- 1 -I, D k = \{I + B k )A k {I-B k ), F a from j25fc). (27c) 

Furthermore, we define A max = max||Efe|| , A~ in = maxUS^" 1 !!, i/j = re lip \ 
-it < cpj < 7T, y = diag(yi, . . . , y p ), 

n 

K = K{ yi ,...,y p ) = Y[ etr(±(Y + B^D^I + BkY- 1 ]-^ , 

k=l 

where the negative sign occurs only with B k , D k from (|27r ). and 

Tadif = F<x(vxj,yj)dipj. 

Theorem 2. With the above notations the functions F from HP are respresentable 
by each of the following three integrals: 

{2n)- p [ KT a dip, (28) 
Jc p 

F a from B k ,D k from {Ffr), \\B k \\ < 1 if v < 2X~l x , max||S fc || < r < 1, 

(j[ etr(-A k )\I + B k \ a >^j {^)" P J c KT a dip, (29) 
F a from Qgfr), B k ,D k from {F7p), \\B k \\ < 1 if v > ±A"^ max||B fc || < r, 

(fl etr(-±A k (I-B k ))\I + B k \ a * J (2w)- p f KT a d v , (30) 
F a from Wfc), B k: D k from (FTjp), v > 0, max\\B k \\ < r < 1. 
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Proof. Because of lemma 1 the assumptions of theorem 1 are satisfied with 
9ja(t) — z f = (1 + 11 ~ 1 tj)~ a an d hj{tj) corresponding to Zj or Uj — v~ 1 tjZj = 1 — zj 
or ujj = Zj — Uj respectively. The functions '§j (t)(hj(t)) n are the L.t. of the functions 
in the second column of (pZ2)) from which type (a) and (c) have the bounds in (j2"3")) , (fM)) . 
satisfying the condition (JT]) for theorem 1. The series J2^=o G a +n{x)y n = G a (x,y) in 
(|25b ) is absolutely convergent for every y 6 C. Thus, all r > max||Bfc|| are admissible 
in (|29p. In (j50)) we have max ||i?fc| < 1 for every v > 0. Hence, theorem 1 together with 
the respresentations of the L.t. in implies fl25J), (HHJ) and P0|) . 

The univariate case of (|29|) provides 

F(x;ai, ...,a n ,a 2 ,.. . ,cr£,tf 2 , . . . ,<5 2 ) = 

with 2u > max cr^ 2 , r = 1, G Q from (|26|) . With p — 1 similar formulas arise from (|28|) 
or (130]!. 

The cdf of a quadratic form \x' Ax with T' AT — diag(Ai, . . . , A n ) > of rank q 
and a Af(fJ,, fj 2 /„)-random vector x is a special case of (pITj) with at = |, er 2 . = A/jCr 2 and 
non-centrality parameters Sf. = hfi^/a 3 , k = 1, . . . , q, /i* = T'/i. 

Some further remarks: In (|29|) also ||Bfe|| > 1 is allowed since every r = ||Y|| > 
max||_Bfc|| is admissible, which entails max ||_BfcF _1 || < 1. 

With i9 = A max /A m i n it follows with special values of v: 

■d — 1 

max||B fc || < —— in (EHD with v = 2(A min + A max )" 1 , 
v + 1 

max||B fc ||<^| in © with v = i(A^ n + A^J, 



but 



maxlLBjUl < ^£ — - in with t> = (A min A rr - ^ 1/2 

1 



Mnax I 



More generally, the scale factor v — w 2 can be replaced by a scale matrix 
W 2 = diag(w 2 , . . . , w 2 ) > 0. Then with T w = W^TW' 1 , T, w = WT,W, A w = \VAW~ 1 
the L.t. ([T]) equals 

\I + -£ w T w \- a eti(-E w T w (I + E W T W )~ 1 A W ). (32) 

Consequently, besides the substitutions wSfc — > WEkW, A& — > VFAfcVF -1 , the matrices 
/ + B fc in theorem 2 must be replaced by W^E fc VK, (VKEfcVTT 1 and 2(7 + WEkW)' 1 
respectively, and the generators F a {vxj ,yj) by F a (w 2 Xj,yj). 

In particular for a single r p (a, S, A)-distribution this more general scaling can be 
used to minimize ||£?|| or for a "natural scaling" i.e. to standardize J + B to a correlation 
matrix. However, ||_B|| < 1 must be taken into account in (I28p . whereas this condition 
is satisfied in (|30|) for every scaling. It was shown in Royen (1991) that natural scaling 
can always be accomplished also in I + B = 2(1 + WSW) -1 by a unique W 2 . 



5. Representations of the T p (a, E, A) distribution function by (p — l)-variate 
integrals 

For a single r p (a, S, A)-cdf it is always possible to perform the integration over a 
single variable tpj within the integrals from theorem 2. 
We use the following functions 

g a (x,y) = e~y J2 G a+n (x) V - = £ G^z)™-, 

n—Q ' n— (^^W 

x,y£ C, G Q+n , G*£l n from and 
g*(x,y)=e^ a ( : r,y). 

For positive half integers a = 1/2 + fc these functions can also be computed by the 
erf-function and a sum of k terms which are essentially given by the modified Bessel 
functions / 7 -_ 1/2 (2(xy) 1 / 2 ), j = 1, . . . ,k, (see e.g. Royen (1995) or (2006)). 

Now let be W 2 — di&g(w 2 , . . . , wT) a general scale matrix, 
Y = diag(yi, ...,y p ),yj = re lH> \ -in < cpj < n, 

, v WYW - 1, (34a) 
B = I u I I = \ (WZW)- 1 - I, (346) 
2{I + WT.W)- 1 - I, (34c) 



h' h 



( v WAZW, (35a) 

D= [ <P d) = \ W " lE " lAW_1 . ( 35& ) 
^ [ 2(7 + WEW^WAW- 1 ^- (J+ WEWO -1 ), (35c) 

Vo = 2/o(2/i, ■ ■ • , 2/p-i) = b' p (Y pp + B pp y\ - b pp (36) 
g = q ( Vl , 8p_0 = (& p (F pp + 73 pp )-\ -1)2? (37) 

and 

K a = K a (y u 2/p _ 1 ) = etr(±(F pp + i? pp ) _1 £' pp )|/ + B^Y^ 1 ]-*, 
where the negative sign is only taken for 73 pp from (134b ). 

Theorem 3. With the above notations the r p (a, S, A)-cdf F(x±, . . . ,x p ; a, E, A) 
is given by each of the following three integrals: 

7^— j— f / G a T ^-) kJyI - J- G a Ux h &p h (38) 

S /rom {3^, D from fffSbj, ||S|| < r < 1, 

eta^-WA^" 1 ) 1 /•,,.. s_ Q / n .. , , 9 



^ Jc^ ( " " o) Qa V " 2/o) " p " p ' T" 



27o 



p-i 



(39) 



^ Q G a (wjXj,yj)d(pj, 
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B from D from \\B\\ < r, 

2 ap etr(—^WAW~ 1 (I — B)) 1 



exp 



\i + wzw\" (27r)f-i 7 Cp _, "Vi-y 0/ (4(|) 

(i - yo)~ a G a ( - Vo w 2 x p , 2q 2 ) # a TT —5— G Q ( w 2 ^, -%-) 

5 /rom {^;, L> from Wc), \\B\\ < r < 1. 

For the proof of theorem 3 the following two lemmas are required. 

Lemma 2. With Y = diag(yi, . . . ,y p ), yj = re lv \ \\B\\ < r, B,D,y a ,q from [5$ , 
A35\) . \3b}) , \31/\ ) the following decomposition is obtained 

etr((y + i?)- 1 J D)|y + J Br Q 

- ctr ((Y pp + B PP )- 1 D PP ) \Y„ + B pp \~ a exp ( —?—) (y P - yo)~ a - 

\yp-yo/ 



Proof. From frequently used formulas for p x p-matrices, (see e.g. complements 
and problems 2.4, 2.7 in chapter lb of Rao (1973)) it follows for 

A nr) b r , 



A = Y + B=\ pp Up 

\b p y p + b pp 

\ A \ = \ A p P \(y P + hp - h 'p A pp b p) = \ y pp + B P P \(y P - yo), 

(4 _1 4- 1 A~ x h h 1 A^ 1 1 A~ x h \ 

PP ^ Vp-yo^PP U P U P A PP y p —J/o PP °P I 

1 ft' 4-1 1 I 

Vp-Vo P PP Vp-Vo / 

and 

trace(A _1 L>) 

= trace ( A p D pp H (A pp b p b' A pp D pp — A pp b p d p ) J H (d pp — b 'A pp d p ) 

V Up - yo J y P - yo 

= trace(^4~ 1 D„„) H — , which implies (|4"T1) . 

y v - 2/0 



Lemma 3. Let be q any number, S r — {y £ — r}, yo an y number with 

?/o | <r, then with F a from h25}) , Q ai Q^ from &33\) , and the negative sign in ±q only for 

■) 



' / ctr (^^\ Fa ( X] y)( y ^ yo )-« y «-l d y 



2iTi Jsr- \y- yo 



S * T^J F <* f rom r < 1, (42a) 

(1 - yo)~ a 0* ((1 - Vo)x, jS^) , F a from (MB), (426) 

exp (^L) (1 - yo)-" g a (i^!, ^) , F Q from r < 1, (42c) 
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Proof. It is sufficient to verify (|42p for the corresponding derivatives f a = -§-F a 
At first, (|32"b) is shown: 

With F a from pSa) and the binomial series for (1 — yo/y)~^ a+n ^ we obtain 



y 3 (e .aw.-) (t 

17 *r \m=0 / \fc=0 



a + n + fe — 1\ / j/o x fc 



y-^dy. 



With 2 = (1 +t) 1 , tt = tz, the last integral has the L.t. 



\m=0 / \fc=0 



Z 

m—n+k 



E " ro ( a+ n t k ~ = - ^r (a+n) - 



Multiplication by (— q) n /n\ and summation over n leads to the L.t. 

z a ( qu \ 1 / t^(1 -ito)*' 

■ ex P ~ i = /-, , /-, r^77 exp 1 



(l-^ ) a V l-uy J (1 + (1 - y )t) a ' \ l + (l-y )t 

and this is the L.t. of #■ £ a f r^-, ^ 

To verify (g2b) we obtain with F a from ([25b): 



r(a + n + fc) f yo^ 



-) l^"" 1 ^ 



r(a + n)fc! \ y 
- V a + (x) T{a + n + k) v k - a x Me xw 

m— n+fc 

= (l - y )- (Q+n) (i - yo)g a+n ((l - yo)x). 

Multiplication by q n /n\ and summation provides (1 — yo) a jf^ G* (jl 
(H2t) can be shown by L.t. in a similar way as (|l2"a ). 



yo)x, j^- o 



Proof of theorem 3. Without loss of generality y p is selected from the variables 
yj = re llfj in Y = diag(y±, . . . , y p ) with any fixed r > \\B\\ . If y p is replaced by a variable 
y with any \y\ then the equation 

\Y + B\ = \Y pp + B pp \(y p + b pp - b' p {Y pp + B pp )- l b p ) = 

has always a unique solution 

y = yo = b p {Y pp + Bpp^bp - bpp 

with 1 3/o | < r since \\B pp \\ < \\B\\. 

Hence, with lemma 2 and lemma 3, theorem 3 is obtained by integration over ip p in 
the integrals of theorem 2 with n = 1. 
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